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Abstract: The next generation of ground-based Imaging Air Cherenkov Telescopes will play an 
important role in indirect dark matter searches. In this article, we consider two particularly promis- 
ing candidate sources for dark matter annihilation signals, the nearby dwarf galaxies Draco and 
Willman 1, and study the prospects of detecting such a signal for the soon-operating MAGIC II 
telescope system as well as for the planned installation of CTA, taking special care of describing 
the experimental features that affect the detectional prospects. For the first time in such studies, 
we fully take into account the effect of internal bremsstrahlung, which has recently been shown to 
considerably enhance, in some cases, the gamma-ray flux in the high energies domain where At- 
mospheric Cherenkov Telescopes operate, thus leading to significantly harder annihilation spectra 
than traditionally considered. While the detection of the spectral features introduced by internal 
bremsstrahlung would constitute a smoking gun signature for dark matter annihilation, we find that 
for most models the overall flux still remains at a level that will be challenging to detect, unless one 
adopts somewhat favorable descriptions of the smooth dark matter distribution in the dwarfs. 
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1. Introduction 

Even though basically nothing is known about its underlying nature, compelling evidence has ac- 
cumulated in recent years that the building block of the observed structures in the universe is a new 
form of non-baryonic, cold dark matter (DM), with the fraction of DM to critical density being 
~ 0.233 ± 0.013 according to the most recent estimates [[T]]. A class of particularly well- 
motivated DM candidates arises in many extensions of the standard model of particle physics in 
the form of weakly interacting massive particles (WIMPs) that, thermally produced in the early 
universe, automatically acquire a relic density matching the observed DM abundance today (for 
reviews, see [|, H g 1]). In the following, we will restrict our analysis to the supersymmetric 
neutralino, which is maybe the best-motivated and the most widely studied DM candidate. 

The search for DM is a vast field where different experimental techniques can provide comple- 
mentary information. At the Large Hadron Collider (LHC), the detection of new supersymmetric 
particles through missing energy signatures, would provide valuable information about the compo- 
sition of the neutralino and its role as a DM candidate (see, e.g. ||^ ^). The neutralino could also 
be observed through nuclear recoils in direct detection experiments. A new generation of cryogenic 
detectors is aiming to improve the already impressive recent limits reported by Angle et al. [j8|] and 
Ahmed et al. |^, thereby probing significant regions of the underlying parameter space. Finally, 
DM can be searched for through indirect detection, i.e. by looking for primary and secondary prod- 
ucts (mainly gamma rays, neutrinos, positrons and anti-protons) of DM annihilation. Our attention 
will be focused on the detection of high-energy gamma rays. 
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In this field, Imaging Air Cherenkov Telescopes (lACTs) play a leading role, exploring the 
possibility of tracing Cherenkov photons from electromagnetic showers produced by gamma rays 
impinging on the Earth's atmosphere. Ground-based gamma ray astronomy is a very active field 
that has matured rapidly in recent years. In the near future, new installations with improved sensi- 
tivity will investigate those regions in the sky where DM signatures can be expected. In particular 
MAGIC II, already in the commissioning phase, will perform observations in the northern hemi- 
sphere, while a new generation of lACTs is currently in the design phase, the Cherenkov Telescope 
Array (CTA), where an array of several tens of telescopes will have improved angular and energy 
resolution, as well as a very low energy threshold and a greatly increased sensitivity. 

Dwarf spheroidal galaxies (dSphs) are the faintest known astrophysical objects believed to 
be dominated by DM [10]. So far, almost all detected dSphs are located in the Local Group, a 
fact that is most likely related to the low luminosity of these objects, ranging from 330 L© to 
3 X 10^ Lq [11, 12, 13, Typically, dSphs are satellite galaxies hosted by the halo of a larger 
and more massive galaxy that, in some cases, strongly influences their evolution: tidal stripping, 
e.g., can deprive a dSph from stars in the outer region (see, e.g., [15, [l^). While dSphs share 
many characteristics with stellar globular clusters — a similar number of stars (from tens to a few 
thousand) and a similar range in luminosity — they usually have a larger size, a fact that implies 
the presence of a DM halo of the order of 10^ - 10^ M©. The resulting mass-to-light ratios are 
typically larger than 10 Mq/L© and can even reach values up to 10^ Mq/Lq, suggesting that these 
objects are particularly interesting targets for indirect DM searches (for previous studies, see e.g., 
[Il71,[l8|,|l|,|^|2l|]). 

Let us mention in passing that the number of detected dSphs is around one order of magnitude 
below what is expected from A/^-body simulations of structure formation in ACDM cosmologies, 
which has become known as the missing satellites problem [22, 23]. While many interpretations 
have been suggested, including the effect of decaying DM (see, e.g., [|^, 25, ^ and references 



therein), the most natural explanation would be that the missing dSphs have so far simply escaped 
detection. In fact, it has been proposed that dSphs below a certain mass cannot efficiently accrete 
baryons, which would explain the intrinsic faintness of these objects [18]. The claimed discrep- 
ancy has recently also been mitigated considerably by the discovery of a bunch of new ultra-faint 
galaxies in the SDSS data 

In this paper, we provide an update on the possible future of the indirect detection of DM with 
lACTs. In particular, we discuss the prospects for the detection of DM in dSphs with MAGIC II 
and make tentative predictions in view of the expected characteristics of CTA. We choose to focus 
on two dSphs, Draco and the recently discovered ultra-faint Willman 1 . While the choice of Draco 
is supported by the large set of available data, Willman 1 will be discussed as a very promising 
target for future detection. 

A key ingredient in estimating the gamma ray flux for lACTs is the computation of the energy 
spectrum for photons produced by neutralino annihilations. In Ref. [27], it was recently shown 
that internal bremsstrahlung (IB) provides an important, hitherto neglected, contribution to the 
spectrum. While the effect of IB is largely model dependent, it appears in general as a pronounced 
"bump" at energies close to the kinematic cutoff at the neutralino mass. The importance of this 
effect is two-fold: first, the flux at high energies, where lACTs are most sensitive, is significantly 
increased; secondly, the introduction of spectral features allows an easier discrimination of a DM 
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source from potential astrophysical sources located in the vicinity, whose spectrum is usually a 
featureless power law. This is the first time that IB is fully included in studying prospects of DM 
detection for lACTs. 

The paper is organised as follows: in Section ^, we will briefly review some of the main prop- 
erties of Draco and Willman 1 and discuss how their DM profiles are constrained by observations. 
The following Section ^ is dedicated to a short description of the lACT technique, with particular 
attention to the features that are relevant for DM observations. In Section ^, we introduce a set of 
five benchmark neutralinos, representatives for the relevant regions of the parameter space, and es- 
timate the corresponding flux based on the technical specifications for MAGIC II and CTA. Results 
are discussed in Section ||, and in Section ^ we present our conclusions. 



2. Phenomenology of Draco and Willman 1 



2.1 Draco 



Draco is located at a distance of 80 ± 7 kpc from the Earth [11] and is one of the best known and 



most often studied dSph. Discovered by Wilson [28], the first estimation of its DM content was 



performed only 30 years later by Aaronson [29] from the analysis of the stellar velocity dispersion. 
Nowadays, a large set of data is available for Draco [30,^,32,33,^]. Even though the formation 
of the DM halo cannot be traced back with great precision, it can be inferred that the system is 
composed of very old, low-metallicity stars, without any significant sign of star formation in the 
last 2 Gyrs. At its short distance from the MiUcy Way (MW) center, Draco has likely been stripped 
off its outermost stars in the past; today, however, no signs of tidal interactions are observed [^]. 
From a kinematical analysis of a sample of 200 stellar line-of-sight velocities (with radii that range 
from 50 pc to 1 kpc) one can infer the DM profile: The result of the fit (assuming a Navarro- 
Frenk-White profile, hereafter NEW [p5|]) indicates a virial mass of the order of lO'' Mq [||], with 
a corresponding mass-to-light ratio of M/L > 200 Mq/Lq that characterises Draco as highly DM 
dominated. 



Many groups have analysed the available data for Draco and modelised its DM profile [ 19, 2C , 
|36| , |37| , |38| , 39]. We are going to discuss here the two extreme cases of a cusp profile and a core 
profile, which span the range of possible configurations.^ Eor the former, we take the NEW profile. 



Pnfw ('') — 



(r/r,) (1 +rlr,f' 



and for the latter we take a Burkert [ ]41| , |42| ] profile. 



P Burkert (r) 



Ps 



(2.1) 



(2.2) 



(1+r/r,) (l+(r/r,)2)' 
Eor the scale radius and the scale density ps, we use the values summarised in Table |l| (taken 



from Ref. [39]). Note that both profiles produce good fits to the velocity dispersion profiles, down 
to the pc scales where the innermost stars are observed. No direct observational information is 
available for distances even closer to the centre, where the two profiles differ significantly. 



'It should be underlined that core profiles are not compatible with current CDM simulations: recent results predict 
DM halos well-fitted by an Einasto profile without converging to a precise value for the inner slope |po|]. For comparison 
with previous studies, however, we chose to still include core profiles in the analysis. 
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Draco-Burkert Draco-NFW 


Willman 1-NFW 


r, (kpc) 


0.35 0.50 


0.18 


p, (Mo/kpc^) 


3.6x10^ 1.3x10^ 


4.0 X 10^ 


090 n 


0.52° 0.35° 


0.20° 


D (kpc) 


80 ±7 


38 ±7 



Table 1: Scale radius (kpc) and scale density (Mg/kpc^) that appear in the DM density profiles. The NFW 
and the Burkert profiles in the case of Draco represent models N3 and B2, respectively, from Ref. [^. In 



the case of Willman 1, the NFW fit is taken from Ref. |52|. In addition, the semi-aperture of the solid angle 



corresponding to 90% of emission flgo (°) is reported together with the target distance D (kpc). 



Possible DM annihilation signals from Draco have already been searched for in the past: after 
the CACTUS experiment had claimed the detection of an excess of ~ 7000 high-energy photons in 



only 7 hours [143[], almost all lACTs tried to reproduce the result, but the claim was not confirmed 



[^, MAGIC observed Draco for 7.8 hours during 2007 above 140 GeV ||4^]. Within a 
2a" confidence limit, the collaboration reported an upper limit for the integral flux of O" ' (£" > 
140 GeV) < 1.1 X 10~i^ ph cm'^ s'l. 

2.2 Willman 1 

Willman 1 (SDSS J1049-I-5103) is a very peculiar object, located at a distance of 38 ± 7 kpc from 
the Earth in the constellation of Ursa Major. Discovered in 2004 by B. Willman pTl Bsl], using 



data from the Sloan Digital Sky Survey [|49|], it was then further observed with Keck/DEIMOS [|50|] 

and more recently by Siegel et al. [51 ]. With an absolute magnitude of My 2.5 and a half-light 

radius (i.e. the radius of a cylinder, pointing to the Earth, that encloses half of the luminosity of 
the object) of 21 ± 7 pc, it looks very similar to a globular cluster, even if its narrow distribution of 
stellar velocities and the large spread in stellar metallicities suggests that it is indeed the smallest 
dSph ever observed. The object may show evidence for tidal disruption from its tri-axial stellar 



distribution [|48|]. On the other hand, the difference in the stellar luminosity function between the 
central and outermost stars reveals a strong mass segregation. With a luminosity of 855 L©, and a 
mass of the order of 5 x 10^ M© [pO|], Willman 1 could feature a mass-to-light ratio in the range 
500 - 700 Mq/Lq, or even more, making it one of the most DM dominated objects in the Universe 



& 

The small number of stars that belong to this dSph hinders, however, an accurate determination 
of the DM density profile. Following Strigari et al. [^], we parametrise its DM halo with an 
NFW profile, as specified in Table [l], although these parameters are subject to somewhat larger 
uncertainties than in the case of Draco. 



3. Observations with lACTs: MAGIC II and CTA 

When entering into the atmosphere, cosmic gamma rays (as well as the many orders of magnitude 
more frequent background charged cosmic rays) quickly lose energy through interactions with 
the nuclei of atmospheric molecules, dominantly by pair production of electrons and positrons. 
Bremsstrahlung photons radiated by these highly energetic electrons and positrons in turn lead to 
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the production of secondary electron-positron pairs, thus triggering the subsequent development of 
a particle shower. When the electron-positron energy falls below Ec ~ 83 MeV in the atmosphere, 
the dominant mechanism of energy loss becomes ionization and the shower rapidly dies off. This 
takes place at an altitude of 8 - 12 km, depending on the energy of the primary gamma ray. This 
cascade of highly relativistic particles causes a flash of UV-blue Cherenkov light, with the greatest 
emission coming from the shower maximum (i.e. where the number of free electrons and positrons 
is maximal), lasting a few nanoseconds and propagating in a cone with an opening angle of ~ 1°, 
slightly depending on the primary energy. The resulting circle of projected light at 2000 m asl 
(the MAGIC telescope altitude) has a radius of about ~ 120 m. If a telescope is located inside 
this Cherenkov-light pool, the light can be reflected from the collecting mirrors and focused onto 



a multi-pixel recording camera. An image reconstruction algorithm [ ]53[ ] then allows the recovery 
of the energy and direction of the primary particle, and determines whether it was more likely a 
hadron or a photon. In this way, it is possible to reject up to 99% of the background, constituted 
mainly by sub-showers generated by charged cosmic ray particles, by muons and by the night sky 
background light. 



This technique was pioneered by the Whipple telescope [54], with several successors cur- 
rently operating, including MAGIC, HESS, CANGAROO-III and VERITAS^ In this paper, we 
will address observational prospects for the upcoming MAGIC II telescope system and for the fu- 
ture generation of lACTs, focusing on the case of CTA. MAGIC II is a stereoscopic system of 
telescopes, composed of MAGIC and a second telescope currently under commissioning on the 
island of La Palma, which will start operation in 2009. The stereoscopic view of two telescopes 



(pioneered by HEGRA []55P), together with the improved technical characteristics of the second 
detector, will allow a general improvement in the overall performance of the experiment, in par- 
ticular in terms of energy and angular resolution, as well as energy threshold. The performance of 
the MAGIC II array was simulated with Monte Carlo (MC) tools by Carmona et al. [Q. CTA, on 
the other hand, is the result of an effort for a next generation Cherenkov observatory with increased 
capabilities: normally, one single telescope can cover 1.5-2 orders of magnitude in energy range. 
With the combined use of many telescopes of 2-3 different sizes, CTA should be able to extend 
the energy range to almost 4 orders of magnitude, from roughly ~ 30 GeV to ~ 100 TeV. The 
experiment is still in the early design phase and the final layout of the array is thus far from defined 
yet; the performance, therefore, is still subject to changes. For this paper, we will refer mainly to 



the work of Bernloher et al. [ |57| ] and several private communications within the collaboration. The 
CTA prototype construction could start in 2010, at least for some of the main components, and the 
final installation is foreseen in 2012-15. 

The performance of an lACT in terms of its prospects to detect a DM annihilation signal 
can generally be characterised by a small number of basic parameters, which are described in the 
following (see also Table |^ for a summary of the characteristics for MAGIC II and CTA): 

• Energy threshold: The energy threshold of an lACT can take slightly different values accord- 
ing to the definition. Hereafter we consider it to be the peak of the reconstructed MC energy 
distribution (other definitions being analysis threshold, trigger threshold, etc.). This value de- 



tokyo.ac.jp/ and ittp://veritas.sao.arizona.edu; 



"See, respectively, ittp : //w w wmagic . mppmu . mpg. de[ http://www.mpi-hd.mpg.de/htm/HESS, http://icrhp9.icrr.u- 



-5- 



pends mainly on the reflector area of the telescope: A larger mirror area allows, in particular, 
to collect more photons from the showers and thus increases the chance of discrimination 
against the night sky background light, in particular for low energy showers. The use of 
a stereoscopic system also plays an important role because it enhances the gamma/hadron 
iglh) discrimination power which is weaker at low energy. The energy threshold changes 
with the zenith angle of observation, and sources culminating high in the sky are preferred. 
Reaching a low energy threshold is an important feature for DM studies with lACTs, both be- 
cause of the increased number of photons and because of the enhanced possibility to observe 
the spectral cutoff even for low-mass neutralinos. Making use of stereoscopic observations. 



MAGIC II will have an energy threshold of 60-70 GeV pq], with possible extension to even 



lower energies with improved analysis techniques and new trigger systems currently under 



development []58[]. This value will be further lowered to at least 30 GeV for CTA. The tele- 
scope acceptance for gamma rays around 30 GeV starts to decrease rapidly, but a very strong 
gamma ray signal could probably even be detected at energies as low as about 10 GeV. 

Energy resolution: The true energy of the primary gamma ray E' is reconstructed on the basis 
of a comparing analysis between the shower image parameters and MC events. The probabil- 
ity to assign, after the analysis, an energy E to the primary gamma ray can be approximated 
by 

1 / {E-E'f\ 

Typical values for the energy resolution e are of the order of 10-30% for lACTs, depending 
on the energy. The reason for such large uncertainties is the combined effect of many sources 



of uncertainties (for a more detailed discussion, see Ref. []46[]). The energy resolution is an 
important parameter when observing spectral features as bumps and cutoffs that can provide 
clear signatures for a DM signal. MAGIC II will have an energy resolution of 15% above 
300 GeV (up to 20% at 70 GeV); for the CTA, this situation could radically improve. Finally, 
let us note that further systematic errors might hide in the absolute energy calibration; the 
recent MAGIC observation of a clear cutoff in the Crab Nebula spectrum^ , when compared 



to a corresponding future observation by Fermi-LAT |59], may allow for the first robust 
calibration of gamma ray energies []60|]. 

Angular resolution: The reconstruction of the direction of a primary gamma ray is performed 
through image analysis. As a result, a gamma ray coming from a direction ifj' will be recon- 
structed to a direction in the sky with a probability distribution that can be fitted to a 
Gaussian function: 

R ,^ 1 / {^'-^? \ ,,,, 



^The Crab Nebula is a supernova remnant that is conventionally taken as reference source for cross-calibrations in 
gamma ray astronomy due to its very stable and intense flux. 
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The standard deviation ■Qr of the Gaussian corresponds to the telescope angular resolution, 
also called Point Spread Function (PSF). As a consequence, any source will appear somewhat 
blurred. The stereoscopic system exploited in MAGIC II will improve the PSF, allowing 
values as low as 0.05°, while for CTA we expect an even smaller PSF. It is hard to predict an 
exact value given the current lack of knowledge of the CTA design, but a realistic value that 



we use for this study is 0.02° (see also [61 1). For extended sources, as in the case of dSphs, 
the PSF plays an important role in the reconstruction of the DM density profile, as discussed 
in the next section. 

Flux sensitivity: The sensitivity of an lACT is usually defined as the minimum flux for a 5 cr 
detection over the background, after 50 hours of observation time and based on at least 10 
collected photons. For operating experiments, the sensitivity can be computed by using real 
data and following Eq. 17 of [^], while for planned experiments the sensitivity has to be 
estimated on the basis of MC simulations and is therefore subject to larger uncertainties. The 
procedure is as follows: a full data analysis is performed on two samples of MC simulations, 
one for gamma ray events and one for background events (basically protons and helium), 
during which a number of parameters ("cuts") is optimised to maximise the analysis quality 
factor Q = ey/ yfefj, i.e. the ratio between the efficiency for gamma rays and the square 
root of the hadron efficiency ("efficiency" refering here to the ratio between the number of 
events passing the analysis cuts and the number of events at MC, input, level). After the 
optimisation, one can estimate the number of hadrons Nh{> E) above some energy E. Given 
the Poissonian distribution of events, a 5 o" detection is obtained whenever the number of 
gamma rays detected is larger than 5 yfNk(>'^. The integrated sensitivity above E is thus 
given by: 

^..(,^).5^EI)i, (3.3) 

A ■ tSQ 

where A is the MC gamma ray simulation area and tso is the time interval corresponding to 
50 hours. 

A'^;j(> E) and ey, and thus the sensitivity, are usually determined assuming a featureless power 
law spectrum of MC gamma ray events of index -2.6. This corresponds approximately to 
the spectrum of the Crab. For this reason, the sensitivity is often also expressed in terms of 
"Crab" units (C.U). In the case of the benchmark neutralinos under study (see next section), 
the gamma ray spectra are usually harder than that of the Crab and no longer featureless; it 
is therefore natural to ask how much this would change the sensitivity. 

To address this question, let us note that the sensitivity mainly depends on the g/h discrimi- 
nation power. The g/h sepai^ation, however, is very efficient at intermediate and large ener- 
gies, where the shower parameters are firmly distinguishable between hadronic and gamma 
events. At energies below ~ 30 GeV, on the other hand, the differences are more subtle and 
the sensitivity is affected. Hence we expect that the differential sensitivity does not depend 
too strongly on the spectrum of the source, unless in the case of rather low energies. An 
exact treatment of this effect would require dedicated studies with MC simulations, which is 
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MAGIC 


MAGIC 11 


CTA* 


Eo (GeV) 


100 


70 


30 


e 


30-20% 


20-10% 


10% 




0.10° 


0.05° 


0.02° 


S(> Eo) (cm-2 s-i) 


5 X 10-11 


1.4 X 10-11 


1.5 X 10-11 



Table 2: Comparison of the peri'ormance of the MAGIC, MAGIC II and CTA* telescopes. Eq (GeV) is the 
energy threshold, e the energy resolution and )?r(°) the angular resolution. The sensitivity S(> Eq) (cm"^ 
s"') is given for a Crab-hke spectrum above the energy threshold. 

*For CTA, the numbers have to be taken as placeholders because the telescope design is not yet fixed. 



beyond the aim of this work. Based on a preUminary MC analysis, however, we generally 
expect that the sensitivity at a given energy will not change by more than a factor of two 
compared to that defined for the Crab. A full analysis of this effect is left for future work. 



4. Flux estimation 

The gamma ray flux from DM annihilation can be factorised into two different contributions: 

The term J(iJ/), also called the astrophysical factor, depends on the DM morphology at the emission 
region. The particle physics factor dO^^/dE, on the other hand, depends on the microscopic 
properties of the DM candidate, in particular its mass and cross section, as well as the annihilation 
modes and branching ratios, that define the gamma ray spectrum. For a given telescope energy 
threshold £0, the integral flux is thus (D(> Eq) - J ■ <D^^(> Eq). 

The abiUty to reconstruct the two factors depends on the experimental performance of the 
TACT; in the following, we will estimate the expected flux for Draco and Willman 1. 

4.1 Astrophysical Factor 

The astrophysical factor depends on the source distance and geometry (as well as the PSF of the 
telescope), but for a given DM profile it does not depend on the particular DM candidate. As a 
consequence, the discussion here remains valid for any generic WIMP candidate. Pointing the 
telescope towards a direction if/ in the sky, and taking into account its finite angular resolution, the 
astrophysical factor is given by: 

m ^^jdnjdA [p\r{A,i/r)) ■ B^Xd)] > (4.2) 

where the angular integration dQ - d(p d(cos 6) extends over a cone centered around i/r, with an 
opening angle a few times the PSF 1?^. The integration over A is along the Une-of-sight, in the 
direction ip, such that r = ^JA^ + - 2DAcos(^'), where D is the distance of the source from the 
Sun and cos(^) = cos(0) cos(t^) - cos(^) sm{6) sm{i^). Defined as above, J{ip) is conventionally 
expressed in units of Mgkpc"^ sr"i or GeV^ cm"^ sr 1. In order to translate it to the dimensionless 
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J 


Draco-Burkert 


Draco-NFW 


Willman 1-NFW 


(GeV^/cm^) 
(Mo/kpc^) 


3.84 X 10" 
8.63 X 10^" 


4.71 X 10" 
1.06 X 10^1 


9.55 X 10" 
2.15 X 10" 



Table 3: Comparison of the integrated astrophysical factor J for Draco and Willman 1, for the profiles 
specified in |l} 

quantity J(t]j) as defined in Ref. one simply has to multiply it by 5.32 x 10"^' GeV^^cm^sr 
(= 2.37 X lO-i^Mo^^kpc^sr). 



Integrating Eq. 12 over the full angular extension of the source gives: 

7^ J da^m^^ J dVp\r), (4.3) 

where the second integral is over the spatial extent of the source. Note that this expression is now 
given in units M| kpc~^ (or GeV^ cm~^) and no longer depends on the telescope PSF. 

Using the DM profile parameters of Table |l|, we now show in Figure |l| the quantity 7((A) for 
Draco and Willman 1 . While the two sources, from this plot, appear similar in terms of their angular 
size ifr, we recall that their spatial extension is quite different: By comparing, e.g., their respective 
scale radii for an NFW profile, one sees that Willman 1 (r, = 0. 18 kpc) is considerably smaller than 
Draco (r^ = 0.50 kpc). In the case of Draco, the cusp and core profiles are almost identical (up to an 
overall normalisation factor of ~ 2) for angular distances larger than about tfr ~ 0.3°, below which 
the cusp profile starts to increase more rapidly. At the center, the two profiles differ by around one 
order of magnitude, the difference increasing with decreasing PSF. Whenever an extended emission 
would be observed, one could thus in principle be able to discriminate between different profiles 
by comparing the flux at different distances from the center. As it becomes obvious from Figure |l|, 
and as already stressed by Sanchez-Conde et al. [|^], the angular resolution of the telescope would 
play an important role in this case. Taking into account the full range of profiles consistent with 



the observational data []52[], we note that the astrophysical factor for dSphs is far better constrained 
than for, e.g., the galactic center, where the uncertainty in the inner part spans several orders of 
magnitude [|64|]. 

Given a telescope PSF of the order of 0.1°, and the expected feebleness of the signal, however, 
the capability of reconstructing the morphology of extended sources is very limited. This is partic- 
ularly true in the case of non-stereo lACTs where the shower direction is reconstructed with less 
precision. Even when making the rather optimistic assumption that a signal could be discriminated 
against the background out to a distance where the annihilation flux is a factor of 3 less than from 
the direction towards the center, the source would appear at a size of only roughly twice the PSF 
for both Draco and Willman 1 , in the case of a cuspy profile. For a core profile, the same measure 
would indicate an apparent extension out to < 0.2°, still well contained in a normal lACT camera 
(~ 3° aperture). As we will see, the expected annihilation fluxes are rather low and we find it 
therefore premature to discuss in depth the possibilities to distinguish between different profiles in 
the way indicated above; rather, we will in the following focus on the total, i.e. integrated, flux. 

Table ^reports the calculation of the integrated astrophysical factor / for the two dSphs studied 
here. For Willman 1, the uncertainty in the DM profile translates into a 95% confidence interval of 
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Figure 1: The J{4i) factor in the case of Draco (upper plot) and Willman 1 (lower plot). Core profiles are 
shown in blue, cusp profiles in red. Thick solid (dashed) lines represent profiles smeared with the MAGIC II 
(CTA) angular resolution. Thin solid lines represent the profiles without smearing. The upper right panel in 
each figure shows a zoom-in of the region close to the center For comparison, we also show the profiles for 
a hypothetical, infinite angular resolution. 



about 8 X 10" GeV^/cm^ < 7 < 4 x 10^^ GeV^/cm^ ||5^]. In the case of Draco, the astrophysical 
factor lies in the range 10" GeV^/cm^ < / < 2 x 10^^ GeVVcm^ A gain, these astrophysical 
uncertainties are rather small when compared to other potential sources of DM annihilation signals 
- but one should keep in mind that our choices of DM profiles are actually quite conservative: 
taking into account the above quoted range of possible values for J that are consistent with current 
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Table 4: Parameters defining the benchmark models and some relevant quantities related to the annihilation 
spectrum, mi/2 and niQ (expressed in GeV) are the uniform masses of gauginos and scalars, respectively. tanyS 
is the ratio between the vacuum expectation values of the two Higgs bosons. Aq (in GeV) is the coefficient 
of the trilinear scalar term and fi is the coefficient of the mass term in the Higgs potential, (in GeV) is the 
neutralino mass, Q^/j^ its relic density and crv|i,=() (expressed in cm-' s"') its annihilation rate today, r is the 
ratio of IB photons over secondary photons (above 0.6 m^) and (Tv\yy, crv\zy are the annihilation rates for the 
J lines. O^'' (expressed in cm^ s ' GeV"^), finally, is defined in Eq. 4.4 and given for MAGIC II (£0 - 70 
GeV) and CTA (£0 - 30 GeV) energy thresholds, respectively. 



observations of velocity dispersions in the dwarfs, one could thus win a factor of up to about 4 (in 
the case of Draco) or 40 (in the case of Willman 1) in the annihilation flux; we will get back to this 
in Section B. 



4.2 Particle Physics factor and Benchmarks 



The particle physics factor in Eq. \. \ is given by: 



r dN'AE') 



dE 



(4.4) 



where crv is the total annihilation rate of the DM particles, the mass of the DM particle, B, 
the branching ratio into channel / and dN'^jdE the corresponding differential number of photons 
per (total) annihilation. The integration over Rf:{E - E') (see Eq. ^Tl] ), takes into account the finite 
energy resolution of the instrument. The total number of photons above some energy Eq of course 
no longer depends on the energy resolution (as long as 1 - E^jM^ » e) and is given by: 



Ny{> Eo) 



2- r 



dN!y{E) 
dE 



dE. 



(4.5) 



Three different contributions to the spectrum can, in general, be distinguished: first of all 
monochromatic y lines, where photons are primary products of neutralino annihilation through 
the reactions yy and xx ~^ [63]. These processes are very model-dependent; while 

providing a striking experimental signature, they are usually subdominant (for a recent analysis. 
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see []27p). As a further contribution, secondary photons are produced in the hadronization and 
further decay of the primary annihilation products, mainly through the decay of neutral pions, 
resulting in a featureless spectrum with a power law like behaviour at small photon energies and 
a rather soft cutoff at m^. These contributions, which always dominate the total spectrum at low 
energies, are highly model-independent and have a spectral shape that is almost indistinguishable 
for the various possible annihilation channels (see, e.g., ||4, 64]). 

For charged annihilation products, a third contribution has to be included; internal bremsstrah- 
lung, where an additional photon appears in the final state. As has been pointed out recently [^], 
these photons generically dominate the total spectrum at energies close to the kinematical cutoff 
at m^. Moreover, they add pronounced signatures to the spectrum that would allow for a clear 
identification of the DM origin of an observed source; viz. a very sharp cutoff at and bump- 
like features at slightly smaller energies. While photons directly radiated from charged final states 
give rise to a rather model-independent contribution [|6^], photons radiated from charged virtual 
particles strongly depend on the details of the underlying DM model. 

As anticipated in the introduction, for what concerns the particle physics content, we will 
restrict ourselves to supersymmetric DM. While the Minimal Supersymmetric Standard Model 
(MSSM) needs more than 100 parameters for its full characterisation, in the following we will 
only consider a constrained version, minimal supergravity (mSUGRA - see, e.g., [56|), where 
gravity is supposed to mediate SUSY breaking and the number of parameters is reduced to 4 plus 
a choice of sign after certain commonly adopted, physically well-motivated assumptions such as 
neglecting CP violating or flavour-changing neutral current interactions and imposing the grand 
unification theory (GUT) condition for the gauge couplings. Mi = 5/3tan^0„,M2 w 0.5 M2. The 
parameters defining an mSUGRA model are: universal masses for gauginos (mi/2) and scalars 
(mo), a common trilinear coupling term Aq in the SUSY breaking part of the Lagrangian, the ratio 
tan yS of the vacuum expectation values of the two Higgs bosons and the sign for the coefficient 
fi of the mass term in the Higgs potential. By solving the renormalisation group equations, these 
parameters, defined at the GUT scale, can be translated into masses and couplings at observable, 
i.e. low, energies.^ 

Even if highly constrained, these models permit a rich phenomenology. From a cosmological 
point of view, one can single out five regions of the underlying parameter space that are particu- 
larly interesting as they correspond to models with a neutralino relic density in accordance with 
the value measured by WMAP: the bulk region at low mo and mi/2, where the mass spectrum con- 
tains light sleptons £ and, as a consequence, the relic density is mainly determined by annihilation 
processes xx ~^ in the early universe (through a f-channel exchange of £); the funnel region 
at intermediate values for mo and mi/2, where m^ « 2m^ and annihilations in the early universe 
are thus enhanced by the presence of the near-resonant pseudo-scalar Higgs boson; the hyperbolic 



"^For the calculation of the low-energy features of these models (i.e. mass spectra etc.), we use version 4.01 of 
DarkSUSY | |67[ ] that relies on the public code Isajet 7.69 [ps]]. One should be aware, however, that these calculations 
are highly sensitive to how the renormalisation group equations are implemented and different codes, or even different 
versions of the same code, may give rather different results (see, e.g., Battaglia et al.,2001;2004b). Typically, the 
qualitative low-energy features of a given model can still easily be reconstructed by allowing for slight shifts in the 
parameter space (defined at high energies). From a practical point of view, this situation therefore does not constitute 
a severe problem as one may always regard a set of low-energy quantities like the mass spectrum, annihilation cross 
section and branching ratios as a valid effective definition of the model. 
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branch or focus point region, where mq » mi/2 and the neutrahno becomes an almost pure Hig- 
gsino, with large annihilation rates into gauge bosons; the stau coannihilation region at large nii/2 
but small mo, where ^ and coannihilations with staus (and usually other sleptons as well) 
are important in determining the relic density; and finally the stop coannihilation region (arising 
when Aq + 0) where w mj. 

In what follows, we choose to work with a set of benchmark models, representative of the 
different regions in parameter space described above. From an experimental point of view, the 
advantage of benchmark models is that they allow a direct comparison between data from different 
experiments (most of the benchmarks that we use have already been extensively studied in other 
contexts) and, in general, a more detailed per case analysis than for, e.g., parameter scans. Our 
particular choice of benchmark models is summarised in Table ^ 

The features of these models that are important in our context are the following: 



/': This model (like the following two) was introduced by Battaglia et al. [p9|], where also its 
phenomenology at colliders was extensively studied. It is a typical example of a model in the 
hulk region. While the annihilation into lepton pairs is strongly suppressed for neutralinos 
with the small velocities they exhibit today (unlike in the early universe), annihilation into 
£^€~y, which does not suffer from helicity suppression [[70|], gives a considerable contribu- 
tion due to the lightness of the sleptons. 

J': This model lies in the coannihilation tail. The sleptons being close to degenerate with 
the neutralino, IB from lepton final states gives even higher enhancement of the flux than in 
the previous case. 

K' : A representative model for the funnel region, where the annihilation dominantly occurs 
through an i^-channel pseudo-scalar Higgs boson. Consequently, the additional emission of a 
photon does not lift the helicity suppression in this case and therefore IB contributions have 
to be subdominant. 



F*: Introduced in Ref. [27] as BM4, this model exhibits a large neutralino mass, as typical in 
the focus point region. In this regime, the chargino is close to degenerate with the neutralino 
(in this case an almost pure Higgsino) and large IB contributions result from charged gauge 
boson final states [^T]]. 



J*: Introduced in Ref. [ ]27[ ] as BM3, this is another example of a neutralino in the co- 



annihilation region, characterised by a particularly large IB contribution. 



We used DarkSUSY, which in its most recent public release 5.0.1 [72], contains a full im- 
plementation of the IB contributions focused on here, to compute the annihilation spectra for the 
benchmark models defined above. Line signals are also taken into account, but they turn out to 
be completely subdominant in the cases studied here (except for model F*). The resulting spectra 
are plotted in Figure ^, both before taking into account the finite energy resolution of the detector 
and for the case of an energy resolution of 10%. The main characteristics of these spectra are also 
summarised in Table ^. 
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10^ 




E[GeV] 



Figure 2: The particle physics factor d<^ jdE, as defined in Eq. XA , for the benchmarks models introduced 
in Section 4.2. The upper panel shows the case of a hypothetical detector with perfect energy resolution, 
and a line width of e ~ v ~ 10"^, while the lower case shows the more realistic example of e = 10%. 
For comparison, we also show the spectrum of the Crab Nebula, taken from Ref. [^] with an arbitrary 
normaUsation. 



5. Results and Discussion 

Combining the astrophysical factor of Table ^ and the particle physics factor from Table |^ we can 
finally make predictions about the expected gamma ray flux above the telescope energy threshold 
Eq. a summary of the results is reported in Table ^ where we also quote the increase in the overall 
flux normalization that would be necessary to meet the required sensitivity for a 5cr detection 
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Draco-NFW 




MAGIC II 


CTA30 


i 


0.75 


(1.9- 10^ 1.3- 10"*, 2900) 


4.7 


(3100,2100,490) 


J' 


0.10 


(1.4- 10^3.2- 10^, 7600) 


0.52 


(2.8-10^,4900, 1200) 


K' 


7.0 


(2000, 2000,470) 


35 


(410,260,61) 


F* 


0.45 


(3.M04, 1.6-10^,3800) 


1.1 


(1.3-10^,2800,670) 


J* 


0.37 


(3.8-104,7400, 1700) 


0.42 


(3.4-10"*, 1200,290) 




Willman 1 




MAGIC II 


CTA30 


r 


1.5 


(9200,6200, 150) 


9.4 


(1500,1000,25) 


J' 


0.21 


(6.9-104,1.6-10^380) 


1.1 


(1.4- 10^ 2400, 58) 


K' 


14 


(990,990,24) 


71 


(200, 130,3) 




0.92 


(1.5-10^^,8100, 190) 


2.2 


(6500, 1400,34) 


r 


0.76 


(1.9- 10^ 3700, 88) 


0.85 


(1.7-10\610, 15) 



Table 5: Expected integrated flux <!>(£'> £0) for our neutralino benchmark models (in units of 10 ph 
cm"^ s"' ), where we used the experimental parameters listed in Table |[ In parentheses, we state the increase 
in the signal that would be needed for a 5cr detection as (B1,B2,B3). Here, Bl is the often cited increase 
that is needed when simply comparing the sensitivity and annihilation fluxes above the telescope energy 
threshold E^). B2 is more realistic in that it gives the corresponding quantity above a certain energy E*, 
depending on the benchmark, where the integrated flux to sensitivity ratio is greatest and B3 is the same as 
B2, yet for the most favourable halo profile consistent with current observations (still not taking into account 
the eflect of substructures, however). See text for further details. 



(referred to as Bl in the table). While it is customary to quote sensitivities and actual fluxes above 
Eq in this kind of analysis, we recall that DM annihilation spectra are rather hard, in particular 
when taking into account possible spectral features at photon energies close to the spectral cutoff 
at the mass of the DM particle. On the other hand, the sensitivity of lACTs is considerably better 
at energies somewhat larger than the telescope energy threshold. We therefore take the projected 
sensitivities for the integrated flux above some energy E* > Eq, using the sensitivity curves as 
provided by Bemloher et al. [57] for CTA and Carmona et al. [56| for MAGIC II and, by comparing 
those to the annihilation spectra, compute the minimal increase in the normalisation that is required 
to see at least part of the DM annihilation spectrum above E* . This is referred to as the quantity 
Bl in Table |; finally, we also state as B3 the corresponding value for the most favourable smooth 
halo profile that is consistent with the observational data (i.e. here we take the upper limit on / as 
discussed in Section H~T|). 



So far, we have only discussed smooth DM distributions. On the other hand, it is well known 
from both theory [73] and numerical N-hoAy simulations [0] that cold DM is expected to cluster 
and thereby to form substructures with masses all the way down to the small-scale cutoff in the 
spectrum of matter density fluctuations, which can be determined to a great accuracy from the 
underlying DM model [75]; if surviving until today, such inhomogeneities in the DM distribution 
would greatly enhance the DM annihilation rate [76]. For the case of typical dSphs, this could 
result in an additional boost of the signal by a factor of 10-100 [21]. Another considerable boost in 
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the annihilation flux could also result from the existence of a hypothetical black hole at the center 



of the dwarfs [|19||. In the most optimistic astrophysical configuration, the required increase stated 
as B3 in Table ^ would thus further be reduced by up to two orders of magnitude. 
Let us now discuss some details of the results: 

• Sources. For Draco, the model-dependent fluxes for the Burkert and NFW profiles are very 
similar, and therefore we presented only the latter in Table |5[ For the astrophysical bench- 
mark profiles introduced in Section |4.l| , detectional prospects for Draco and Willman 1 only 
differ by a factor of around 2, and are obviously not very encouraging. When considering the 
most optimistic astrophysical configurations, adopting the highest observationally allowed 
value for J, things change considerably and Willman 1 becomes an interesting and indeed 
very promising target for DM searches. Allowing for an additional, in fact well-motivated, 
boost due to the presence of DM substructures in the dwarfs, this may give at least CTA the 
chance to see also Draco in some cases. 

• Telescopes. Depending on the DM model, the ability of CTA to detect gamma rays from 
DM annihilation is a factor of 6 - 8 better than for MAGIC II. Focusing on Willman 1 , and 
assuming very favourable astrophysical conditions, CTA would in principle be able to see 
all the benchmark models considered here, while MAGIC II should be able to see at least 
some of them. We recall that the flux enhancements needed for a 5cr detection, as states in 
Table ||, are calculated with respect to an observation time of tots - 50 hrs and scale like t~ll.^- 
For prolonged observation times, one could thus win a factor of a few for both telescopes. 
Furthermore, as the CTA parameters are still quite preliminary, an additional factor of 2 in 
the sensitivity of the operating instrument seems quite feasible. 

• Benchmark models. The best prospects for detection are found for the neutralino in the funnel 
region (model K'), the reason simply being a rather large annihilation rate. The second-best 
prospects are found for model /* in the coannihilation region. Recalling from Table ^ that 
/* is actually the model with the smallest annihilation rate, this may come as some surprise 
and nicely illustrates the importance of including IB contributions when estimating the flux 
from DM annihilation. The model F* is yet another example with rather pronounced IB 
contributions; a mass of almost 2 TeV, however, efficiently suppresses the annihilation flux 
(in this case, the required boost actually depends significantly on the details of the — so far 
not sufficiently well known — integrated sensitivity of CTA for TeV photons and could thus 
eventually be significantly improved). 

When compared to previous work, we thus arrive at considerably more optimistic conclu- 
sions than what has been reached by Sanchez-Conde et al. [ pO[ ] for the analysis of present-day 
gamma ray telescopes - not the least due to our fully taking into account all the contributions 
to the expected annihilation spectrum. On the other hand, we find the conclusions of Strigari 



et al. 05 2| ] overly optimistic, a fact that we trace back to the very large particle flux factor of 
= 5x 10^^^cm^s"^GeV~^ that the authors assumed as a fiducial value (this should be compared 
to Table ^ and the corresponding values for our benchmark models, which represent typical neu- 
tralino DM candidates). While it may indeed be possible to find DM models with higher gamma 
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ray yields than considered here, we recall that there exist rather tight general bounds on the allowed 
annihilation cross section and the number of high-energy photons that are produced |77]. Finally, 
we would like to remark that it is certainly promising to combine the DM searches in gamma rays, 
as described here, with observations of the same dSphs at other wavelengths (see also [|19|, [78|]). 



6. Conclusions 



In this article, we have computed the prospects of detecting gamma rays from annihilating DM 
in two nearby dwarf galaxies, Draco and Willman 1, for the upcoming ground-based Imaging 
Atmospheric Cherenkov Telescopes MAGIC II and the CTA telescope array (the latter still being 
in the early design phase). We have focused our analysis on a set of five benchmark models, 
representatives for the parameter space of neutralino DM in the mSUGRA framework, and paid 
special attention to describing those telescope features that are most relevant in this context. For 
the first time in this kind of analysis, we have fully taken into account the contributions from 
radiative corrections that were recently reported by Bringmann et al. [p7[]. As it turned out, in 
fact, taking realistic DM spectra has an important impact on the analysis and, although common 
practice, it can be a rather bad approximation to simply assume a featureless DM spectrum like 
from bb fragmentation and/or to only focus on the total flux above a given energy threshold Eq 
in these kind of studies. The basic underlying reason for this is that realistic DM annihilation 
spectra show a harder energy dependence than the sensitivity of lACTs. Once detected, clear 
spectral features would, of course, have the additional advantage of providing a rather fool-proof 
way of discriminating DM spectra against astrophysical background sources - which is even more 
important in view of the still rather large astrophysical uncertainties involved. 

Although these effects do provide a considerable enhancement of the detectional prospects, 
the expected flux from dSphs remains at a level that, for conservative scenarios, will be challeng- 
ing to detect with the next generation of lACTs. This, rather than the angular resolution of these 
instruments, is the reason why the potential of lACTs to discriminate between different DM pro- 
files in dSphs is limited even in the case of the detection of an annihilation signal; the eventual 
disentanglement between cored and cuspy profiles is probably more promising to perform at other 
wavelengths. 

On the other hand, if one adopts the most optimistic astrophysical configurations that are 
compatible with current observational data of Willman 1, i.e. a favourable DM profile and an 
(9(10 - 100) flux enhancement due to the existence of substructures, all of our benchmark models 
approach the reach of at least the CTA which, for the models studied here, is a factor of 6 - 8 
more sensitive to the annihilation signal than MAGIC II (this is, of course, independent of the 
source). The most promising case of our analysis turns out to be a neutralino in the funnel region, 
characterised by no sizeable IB contributions to its spectrum but a rather large annihilation rate; the 
second best case is a neutralino from the coannihilation region, making up for its small annihilation 
rate with enormously large radiative corrections. 

Having demonstrated that the prospects of indirect DM detection through gamma rays do de- 
pend on the details of the annihilation spectrum, and thus the underlying particle nature, it would be 
interesting to perform similar analyses also for other targets of potential DM annihilation. Another 
further direction of extending the present analysis would be to perform a full scan over the param- 
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eter space of viable models. Finally, we have stressed that the very concept of sensitivity of an 
lACT depends on the spectrum that is observed; in the context of DM searches, this is particularly 
important as DM annihilation spectra can significantly deviate from the usually assumed Crab-like 
spectrum. While we have provided a first estimate of how to proceed in such a case, it would be 
warranting to perform a dedicated analysis, using the full power of state-of-the-art MC tools, in 
order to accurately determine the importance of this effect. 

Finally, we would like to mention that even in the case of negative detection, lACTs could in 
principle put interesting upper limits on the flux which in turn would translate into constraints on 
the combined space of astrophysical and particle physics parameters. Though much smaller than 
for other sources like, e.g., the galactic center, the main uncertainty in this case lies in the overall 
scale of the flux as determined by the details of the DM distribution. This, unfortunately, will 
therefore greatly obstacle any stringent constraint from null searches on the particle physics nature 
of DM for quite some time ahead. 

To conclude, nearby dwarf galaxies - and in particular Willman 1 - are very interesting and 
promising targets for DM searches with the next generation of lACTs. An excellent performance 
of these experiments, in particular in terms of the sensitivity at energies slightly below the DM 
particle mass, will be paramount in such searches. In fact, given the low level of fluxes involved, 
a factor of 2 in sensitivity might decide whether a signal will be seen or not. Complementary to 
such demanding requirements on the experiments, the above discussion should also have made 
clear that it will be very important to collect more astrophysical data and to improve the theoretical 
understanding of how DM is distributed in order to reduce the still unpleasantly large astrophysical 
uncertainties involved. 
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Note added: After the completion of this work, Segue 1 has been presented as yet another ultra-faint Milky Way 
dwarf satellite galaxy, with an expected DM induced gamma-ray flux almost twice that of Willman 1 [[l4[]. Our con- 
clusions can be easily applied to this interesting newly discovered target by simply scaling the boost factors reported in 
Table ^ accordingly. We also mention that the MAGIC collaboration recently reported an upper limit on the observation 
of Willman 1 of the order of lO^'-ph cm^^s"' above 100 GeV using the analysis method described in [F9|. 
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